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Abstract 

We study the generation of two-mode entanglement in a two-component Bose-Einstein conden- 
sate trapped in a double-well potential. By applying the Holstein-Primakoff transformation, we 
show that the problem is exactly solvable as long as the number of excitations due to atom-atom 
interactions remains low. In particular, the condensate constitutes a symmetric Gaussian system, 
thereby enabling its entanglement of formation to be measured directly by the fluctuations in the 
quadratures of the two constituent components [Giedke et al., Phys. Rev. Lett. 91, 107901 (2003)]. 
We discover that significant two-mode squeezing occurs in the condensate if the interspecies inter- 
action is sufficiently strong, which leads to strong entanglement between the two components. 
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I. INTRODUCTION 



Soon after the experimental realization of Bose- Einstein condensates (BECs), rich phys- 

nnn 

ical phenomena have been observed and predicted as well I2, |3|| . In particular, there has 
been a surge of interest in the quantum tunneling dynamics of BECs trapped in multiple 
wells, and much attention has been focused on Josephson effect in such systems 0, . On 
the other hand, it was shown that multi-particle entanglement can be generated in BECs via 
the coherent interactions between the atoms P|, and thereafter BECs have played a promi- 
nent role in the field of quantum information. For example, both multi-particle entanglement 
(i.e. spin squeezing) and two-mode entanglement can be generated in a spin-1 condensate 
with three hyperfine sublevels 0, 3]- In this case, two-mode entanglement, describing the 
inseparability between the two modes (respectively with spin projection m = ±1), can be 
used as quantum information protocols to facilitate quantum teleportation of continuous 
variables 0, 0] . 

On the other hand, in addition to BECs with spin and 1, BECs with two internal 
degrees of freedom, e.g. the \F = l,m = —1) and \F = 2,m = 2) sublevels of ^^Rb, are 
also achievable experimentally and are often termed as two-component condensates 10|, , 
which give rise to novel features such as phase separation 12 1 and the cancellation of mean 
field energy shift [13|| . Such condensates can be viewed as collections of interacting spin-half 
particles and can consequently display multi-particle entanglement through the emergence 
of spin squeezing Q|. Likewise, twjj different kinds of atoms (e.g. ^^K and ^^Rb) can also 
form stable two-component BECs 15]. More interestingly, the interspecies interaction of 
these two-component BECs can be varied by the application of magnetic control thus 
hastening various experimental and theoretical investigations in this field. 

Recently, there have been several discussions on the dynamics of a two-component con- 



densate trapped in a double- well potential [17|, |18[ . The behaviors of such systems, which can 
be realized experimentally with the current technology, are arguably much richer than those 
of single-com pon ent condensates because of the interspecies interaction. For example, in a 
recent paper [isf we studied the tunneling dynamics of a two-component condensate whose 
two components are initially separated by the potential barrier between a double-well. We 
found that in the strong scattering regime atoms in the two components can tunnel through 
the barrier in a correlated manner '^]. As both the scattering strength and the tunneling 
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strength between the wells can be tuned independently with various experimental techniques 
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19|], it is expected that such phenomena will become observable in the near future. 
In this paper, we consider a two-component condensate trapped in a double-well poten- 
tial. Both components of the condensate are initially prepared in the ground state of the 
double well and hence are separable at time t = 0. The objective of the present paper 
is to study and quantify the generation of two-mode entanglement in such a condensate. 
Two-mode entanglement is commonly attributed to the inseparability of the density matrix 
describing two systems, which can be found in various experimental situations and is ver^ 
useful in the applications of quantum measurement and quantum information 



1 and is very 
21, Q. 



In fact, entanglement between two ensembles of atoms with spin has recently been achieved 



22| |. and the concept of two- mode entan- 



experimentally via interaction with polarized light 
glement has been generalized to describe such spin systems [6|]. Since each component of a 
two-component condensate in a double-well can be described as a collection of spin-half par- 
ticles, with the two spatially localized modes in the two wells playing respectively the roles of 
spin-up and spin-down states, the condensate is in fact equivalent to two interacting gigantic 
spins that can be described in terms of continuous variables Q. While it has been shown 
that the intra-species interaction is able to create spin-squeezing for a single-component 
condensate trapped a double- well potential in this paper we will demonstrate that the 
inter-species interaction is responsible for the generation of two- mode entanglement. It is 
the intriguing interplay of the inter-species and the intra-species interactions that sparks our 
investigation in such systems. By applying the Holstein-Primakoff transformation (HPT) 
[2^, which reduces this two-spin system to two coupled oscillators, we show that our sys- 
tem is exactly solvable as long as the number of excitations due to atom-atom interactions 
remains low. More remarkably, the condensate in fact forms a symmetric Gaussian system. 
Therefore, we can directly evaluate the entanglement of formation (i.e. the von-Neumann 
entropy for a pure state) of the system from the fluctuations in the quadratures of the two 



constituent components 



21 



24l |. In accordance with the schemes proposed recently 



by Duan et al. [20, |2jj and Giedke et al. [2J], we define a two-mode entanglement pa- 
rameter that measures the fluctuations in the quadratures and in turn analyze the degree 
of entanglement between the two components of the condensate in the present paper. Our 
discovery is that strong interspecies interaction can lead to signiflcant two-mode squeezing 
in the condensate and hence strong entanglement between the two components. 
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The structure of our paper is as follows. In section II, we introduce the Hamiltonian of 
our system and classify it into symmetric and asymmetric cases according to the properties 
of the constituent condensates. In section III, we consider specifically a system with a large 



number of atoms. An effective Hamiltonian, which is exactly solvable, is then derived 
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the HPT. In section IV, we introduce the two-mode entanglement parameter 
to describe and quantify the entanglement in our system. In section V, we study in detail the 
two-mode entanglement parameter for several typical cases. Finally, we discuss the physical 
meaning of the two-mode entanglement parameter in section VI and consider generalization 
our approach to mixed states. 

II. TUNNELING TWO-COMPONENT BEC 

We first consider the tunneling dynamics of a two-component BEC trapped in a symmetric 
double-well potential. The total number of atoms in components A and B of the condensate 
are Na and Ni, respectively. We will further assume that the interaction between the atoms is 
sufficiently weak and adopt the two-mode approximation to describe the tunneling process. 
Under such approximation, the condensate dwelling in each potential minimum is adequately 



described by a single localized mode function j25 
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. Considering the effect of quantum 



tunneling and the conservation of the particle number of each component, we obtain the 
Hamiltonian of the system: 



+— 
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[ala^f + [a^anf] + ^ [{b{hf + [b^^bnf] . (1) 



Here d] {a-) and ft] ib^) are respectively the creation (annihilation) operators of compo- 
nents A and B residing in the j-th well, j = L, R. Since there are two spatial modes (the 
L and R modes) available for each component, the Hamiltonian above in fact consists of 
four bosonic operators. Besides, the parameters Qai^b), i^a{f^b) and k are the tunneling, in- 
traspecies interaction strength of component A{B) and the interspecies interaction strength 
respectively. 

For the convenience of the subsequent discussion of two-mode entanglement, it is instruc- 
tive to represent this Hamiltonian in terms of angular momentum operators by following 
through the standard Schwinger oscillator model to construct a set of spin operators for 
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each component 28|: 
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relations. In the following we will denote the eigenstates of J^ = J^^ + J^y + J^^ and 

Jaz with \ja,ma) SUch that Jl\ja,ma) = ja{ja + l)\ja,ma), and Jaz\ja,ma) = ma\ja,ma), 

where ja = Na/2. In terms of 3a and J;,, the Hamiltonian can be rewritten as 

H = VtaJaz + ^bJbz + l^aJ ax "I" ^b-Jbx ~^ '^f^Jax-Jbx- (3) 

In the absence of atom-atom interactions, the ground state of the system is obviously 
given by the product state \ja, —ja) \ jb, —jb) and the two components are not entangled. In 
the following discussion we consider how the atom-atom interactions affect the evolution of 
the initial state 

\-^it = 0)) = \ja,-Ja)\jb,-3b), (4) 

and show that the two components will get entangled through the inter- and intra-species 
interactions. 

To facilitate later discussion on the phenomenon of entanglement, we further classify two- 
component condensates according to the symmetry properties of the two components consti- 
tuting the condensate. In a symmetric two-component BEC, the parameters of component A 
and component B are equal to one another, namely Qa = ^b and Ha = i^b- These conditions 
hold approximately for condensates consisting of the hyperfine states \F = l,m = —1) and 
\F = 2,m = 1) of ^''Rb llj. The two components of such condensates have essentially same 
masses and magnetic moments and hence Q n = ^b- Their intra-species scattering lengths 
are quite close and, in addition, Ka ^ ki, ^ k [10| . In a quasi-identical two-component BEC 
where Qa = ^b = ^ and Ha = i^b = i^, the Hamiltonian Q reduces to: 

Hi = nj, + kJI (5) 

where J = Jq + Jf, is the total angular momentum of the system. Despite that the Hamil- 
ionian of such a two-component condensate is identical to that of a single- component one 
1^ . the distinguishability of the two species entails the study of two-mode entanglement. 



Meanwhile, for asymmetric two-component BECs relevant physical parameters of the 
two components are generally different. For example, it has recently been observed in the 
experiment that the condensates of potassium and rubidium (Rb-K), which Jiave different 
scattering lengths and masses, can form stable two-component condensates [15]. Therefore, 
it is deemed appropriate to develop a generic analytical scheme to study such condensates. 

n 

In the following discussion, we will make use of the HPT |23] to carry out a thorough 
analytical investigation on the entanglement between the two constituent components. 



III. BOSONIC OPERATOR APPROXIMATION 



In this section, we consider the evolution of a condensate with large numbers of atoms 
and sufficiently weak scattering strengths, namely Na(b) ^ 1 and Qa{b)^f^a, i^b, f^- As the 
initial state, given by (jlj), is the ground state of a non-interacting condensate and the 
scattering strengths are weak, only the low-lying states will be excited in the evolution 
and the coherence of tunneling process can be maintained. The current situation is in 
contradistinction to our previous study that discovered correlated tunneling of the two 
components in the strong scattering regime. However, we will show that the inter-species 
interaction does lead to nontrivial entanglement of the two components. 

To proceed, we apply the HPT to map angular momentum operators into bosonic oper- 
ators I23I . I29I and show that under the HPT our system is in fact equivalent to two coupled 
harmonic oscillators. In HPT, the angular momentum operators 

and Jaz are expressed in terms of bosonic operators a, d^: 



Ja+ = a^j2ja - a'fa, (7) 



Ja~ = \^\/2ja- a^aj a, (8) 
Jaz = (a^'a-ja)- (9) 

Here d^ and d are standard bosonic operators satisfying [d, d^] = 1. Hence, the Hamiltonian 
^ can be written as 

2" 



K,aJa I L d^d / a^a 



-a 
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a'a 



Since Vtao,) ^ ^a, ^^6, it is arguable that 



< 1, 



(10) 



(11) 



leading to an approximate effective Hamiltonian, 

i/eff = ^aa^a + Vtjj'^h + ^ [/t„ja(a^ + a)^ + ^^b3b{y^ + ^)^ + 2/ty^(a^ + a) (6^ + 6)] . (12) 

This Hamiltonian is analogous to that of two coupled oscillators and completely captures 
the essence of the dynamics of the two interacting BECs. Correspondingly, the initial state 
dl)) is given by the vacuum state |0a, Of,) of the two decoupled oscillators described by the 
first two terms in i^ofr, where |na,nfe) represents the Fock state of the oscillators. 

To study the two-mode entanglement in the condensate, it is advantageous to make use 
of the position and the momentum operators: 



Qb 



V2 

-^(6t + 6), p, = ,-L(5t_S) 



(13) 
(14) 



and to rewrite the effective Hamiltonian as 



HeS ="-y(€+ Pi) + y(^6 + Pi) + l^aiaql + ^^bibql + '^H^ jajbqaQb- (15) 

It is then straightforward to solve the resulting equations of motion of qa{t) , Pa(t) , qb(t) and 
Pb{t). For convenience, we express the solutions in matrix form, which reads 



X{t) = {qait),Pa{t),qb{t),Pbit)f 
= U{t)X{t = 0). 



(16) 



Here U{t) is a real 4x4 matrix representing the evolution operator and can be written as 



U 



C El 
E2 D 

with the 2x2 matrices C, D and E being explicitly given by 



v 



c 



1 



/i2 COSCJit — jli COSCo'2t fia(/^2 siu CJit/cJi — jli SmUJ2t/uJ2) 

(/i2 — /ii) \ — (a;]^yU2 sincuit — u;2/ii sinu;2^)/^a /i2 cos cjit — cos ^72^ 



D 



1 



El = 



and 



(/i2 - /il) 



— /il COStUit + yU2 COSCo'2t — fife (/il sin CJit/coii — yU2 sin Co'2i/'^2) 

(/iicui sin cuit — fi2^2 sin ciJ2t) / fib —j^i cos c^it + /i2 cos uj2t 

— (coscuit — cosct;2t) — fib(sinco'it/c<Ji — sinco'2^/^2) 
uji sincjit — UJ2 sinci;2t)/fia — fife(cosct;it — co^i uj2t) / Via 



— r2fe(cosu;it — cosuj2t) /Via — fife(sinu;it/ct;i — sinci;2i/^2) 
{uJi sin cjit — UJ2 sin ci;2t) /fia — (cos Uit — cos c<j2t) 




where 



/il(2) 



1(2) 



and the normal mode frequencies of the coupled oscillation are 



(17) 



± 



1/2 1/2 



(18) 



This approximate solution, which is based on HPT, is valid as long as the condition (jlip 
holds. However, if the normal frequency is complex, the system will become unstable. We 
can therefore determine the stability condition from (fTH|l : 



fife 



— + 2Ka]\— + \ = Kc 



(19) 



the HPT fails to yield a self-consistent solution for systems violating the inequality. It is 
interesting to note that in the limit where ja,jb ^ oo, the condition for stability reduces to: 

(20) 



< a/ l^al^b 

which is a well known result for BECs in extended space, and violation of (12011 will lead to 

n n 

the phase separation of two component BECs |12l . I30l | . It is also worthwhile to note that 
the stability criterion ()19|) depends on the numbers of atoms in the double- well and similar 
dependence has previously been found for two-component BECs in a single well [s^. We 
will, however, assume the stability criterion (fT^ is satisfied throughout the present study 
and obtain analytically the two-mode entanglement parameter for the condensate, which 
will be defined in the following section. We will see that in addition to yielding the analytic 
solution to the tunneling dynamics, the HPT performed here also facilitates our study on 
two- mode entanglement. 
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IV. THEORY OF TWO-MODE ENTANGLEMENT 



Entanglement between two systems that are described in terms of continuous variab 
is usually indicated by an inequality in its Einstein-Podolsky-Rosen (EPR) uncertainty 



es 



^{([A(ga + g.)f) + ([A(p,-p,)f)} < 1, (21) 

where [qm,Pn] = i^mn for m,n = a,b. qa(b) and Paip) ai'e respectively the position and 
momentum operators (or any pairs of quadratures) of system a(6), and the above inequality 
simply implies that the positions (momenta) of the particles are strongly anti-correlated 
(correlated). In general, condition ()21|) is only a sufficient condition for entanglement and 
does not imply separability of the two systems even if it is violated 0, . However, it has 
recently been shown that a necessary and sufficient condition for entanglement, which is 
analogous to (j2II), can be established if the combined system is a Gaussian one in the sense 
that if its Wigner characteristic function, defined by: 

X^'^Aa, A,) = tr[pexp(A,a - Ka) + \h - \lb^)] , (22) 



is a Gaussian function of Xa and Af, j20, |2J|. Without loss of generality, one can assume 
that the expectation values of all quadratures vanish and hence the Wigner characteristic 
function of a Gaussian system is expressible as: 

1 



X^'"KK,Xb) = exp 



-A^MA 
2 



(23) 



where M is a 4x4 real symmetric matrix and the matrix A is defined by A = (A^, X^, A^, X^)'^. 
As the characteristic function can also be written as: 

X^'^HXaAb) = tr[pexp(zV2A^X)] , (24) 

it is obvious that the matrix elements of M are the correlation functions of the quadrature 
variables X = {qa^VaAb^Vb)^ ■ In fact, Mij = {{XiXj + XjXi)) and therefore M is termed 
the covariance matrix. 

As the amount of entanglement between the two systems is unaffected by local unitary 
operations, say local rotations in the q-p plane and local squeezing operations, the matrix 
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M can be transformed into a standard form by several local operations 

ci 0^ 



0: 



n2 C2 
ci mi 

C2 777,2 



where ni,n2,mi,m2 are positive numbers, and 

774 — 1 772 — 1 



777i — 1 7712 — 1 ' 

|Ci|-|c2| = \/ (771 - l)(77li - 1) - ^^(772 - l)(77l2 - 1). 



(25) 
(26) 



It has recently been shown by Duan et al. that a Gaussian system is entangled if and only 
if the following inequality is satisfied 0, Q| '■ 



alijii + 772) - 2(|ci| + IC2I) + (777i + r772)/ao < 2al + 2/a, 



0; 



(27) 



where 



/r77i — 1 



/r772 - 1 



" V "-1 - 1 V ^2 - 1 

We therefore accordingly construct a two-mode entanglement parameter 



(28) 



2G 



2J: 



6 



ao(77i + 772) - 2(|ci| + IC2I) + (ttii + 7n2)/c 



(29) 



2al + 21 al 

Physically speaking, > is merely a suitably weighted EPR-type uncertainty in the 
squeezed quadratures of the system. The sufficient and condition for entanglement men- 
tioned above can then be expressed in terms of an inequality involving this parameter: 



(30) 



If, in addition, the system is a symmetric one such that 771 = 772 = ttii = 7772, Giedke 



et al. j2^ have recently shown that the parameter can also yield the entanglement of 
formation (EOF) of the system, Ep, - 



Epiit) = c+(ei) log[c+(6)] - C-(6) log[c-(6)] if < 6 < 1, 
= otherwise, 



(31) 
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where the functions c±{z) are defined by c±{z) = ± (z)^/^]^/4. It is noteworthy that 

EOF is a proper measure of the degree of entanglement between two systems and is equal 
to the von-Neumann entropy if the compound system remains in a pure state. However, 
unlike the von-Neumann entropy, EOF still works for mixed states. In fact, Josse et al. have 
recently measured the EOF of a pair of non-separable light beams with this scheme js^. 
In the following section, we shall make use of the two-mode entanglement parameter to 
study how the two components of a BEG condensate trapped in double-well are entangled. 



V. TWO-MODE ENTANGLEMENT IN BEGS 



To study two-mode entanglement in BECs, we first show that a two-component BEG 
indeed forms a Gaussian system. Since 



X^'^KKAb.t) = tr{pexp[zy2A^X(t)]} 

= trjpexp iV2{u{tflCfx{t = Qi) 



(32) 



x'^^\^a,^b,t) can be obtained from x^^\^a,^b,t = 0) and we therefore consider the 
characteristic function at t = 0. For the initial state \0a,0b), it is readily shown that 
X^'^K'^a, Afc, t = 0) is a Gaussian function and M{t = 0) is an identity matrix. x^^K^a, t) 
at other times can simply be obtained by the replacement A U{t)A. Gonsequently, we 
show that x^'"KK,h,t) is a Gaussian function of A and: 

1 



expi-- 



U' {t)A M{0)U' {t)A 



(33) 



which directly yields the covariance matrix M{t): 

M{t) = U{t)M{0)U^{t) 
= U(t)U^it). 



(34) 



Hence, the condensate constitutes a Gaussian system to which the two-mode entanglement 
parameter applies. 

After applying several local unitary transformations to M[t), we obtain the matrix Ms[t) 
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of our system, which reads: 

M,(t) 



^ n c 0^ 
On -c 
c n 

y — c n J 



Moreover, the matrix elements of the matrix Mg are expressible in terms of the variances of 
the physical quantities as follows: 

n = 2{(g^(t))(p^(t))-[Re(g„(tK(t))]2}'^' 

= 2{{qm){pl{t)) - [MQbit)p,it))rf\ (35) 
c = 2 [Re{{qa{t)pt{t)){qb{t)pait)) - {qait)qtit)){p,{t)ptmf' , (36) 

and the variances are given in the appendix for reference. It is then obvious that our system 
is a symmetric one with ni = n2 = mi = 1712 = n and hence = 1. Therefore, one can use 
the corresponding two-mode entanglement parameter, which is given by 

it = n-c, (37) 

to evaluate the EOF of the system [2^. 

Now we are ready to investigate the entanglement in our system with the help of 
and Ep{it). We first consider a quasi-identical two-component condensate (e.g. a Rb-Rb 
condensate), where k = Ka = Hb and Via = VL^ = Vt. The time evolution of for such a 
condensate is shown in Fig. ^ where Na = = 400, Qa = = 50, = = 1 and 
K = 0.50 (Fig.IHa)); k = 1.00 (Fig. Cfb)); k = 1.12 (Fig. C^c)). Hereafter we adopt suitable 
units in which Kb = I for purpose of convenience. The solid line and empty circles are 
respectively results obtained from the HPT and numerical diagonalization of the original 
Hamiltonian, showing that the HPT indeed yields a good approximation in this regime. 
We also show the time evolution of the EOF, which equals the von-Neumann entropy in 
this case, by the dashed line in Fig. ^ It is obvious that our system is able to generate a 
substantial amount of two-mode entanglement for most of the time. Besides, one can see 
that there is a strong anti-correlation between these two curves, which can be understood 
as Epiit) is a monotonically decreasing function of C,t for < < 1. 

It is remarkable that the degree of entanglement depends crucially on the inter-species 
interaction. Of course, it is obvious that the EOF or the entropy is zero when the interspecies 

12 



interaction k vanishes. As shown in Fig. CJa), (b) and (c), the entanglement parameter 
decreases while the EOF increases with increasing k. Therefore, one can achieve optimal 
squeezing by properly controlling the interaction parameters. Besides, it is worthy of remark 
that in Fig. ^c) > k > Kg. Therefore, the system is still a stable one and the HPT 
remains valid. Physically speaking, the availability of the two spatially separated modes in 
fact stabilizes the condensate despite that k > js^. 

We now switch our attention to asymmetric BECs consisting of two components with 
different physical characteristics. In fact, stable BECs of rubidium and potassium have 
recently been achieved in experiments and it is also possible to control the strength of 
interspecies interaction between the two components with a magnetic field lift]. It is therefore 
deemed appropriate to investigate how the entanglement in such condensates changes with 
the interspecies interaction strength, k. In the following, we assume that the tunneling 
strengths and the intraspecies interaction strengths of the two species are in the ratios of 
1:1.45 and 1:1.33 respectively, which are reasonable estimates of experimental data for a 
Rb-K condensate [2^. 

The two-mode entanglement parameters for three asymmetric cases with k = 
0.5,0.875,0.965 are respectively shown in Fig. Efa), (b) and (c). Considering the inter- 
species interaction n as an adjustable parameter, we find that a smaller (i.e., higher 
entanglement) can be obtained as the system becomes closer to the point of stability limit 
given by (|T9|) . As shown in Fig. El^c), if k is increased and approaches the stability limit 
given by (fT!I|l from below, the two-mode entanglement parameter (the EOF) can attain much 
smaller (greater) values. Thus, the significance of the strength of interspecies interaction in 
two-mode entanglement generation is clearly demonstrated. 

From the results illustrated in Figs. 1 and 2 it is manifest that a substantial increase 
of entanglement can be achieved in the vicinity of the stability limit ()19|). In fact, this 
novel feature can be understood heuristically as follows. In general, position (momentum) 
squeezed states of a harmonic oscillator can be produced by strengthening (weakening) its 
spring constant 33j . The two-component condensates can be viewed as a coupled oscillators 
and at the critical point of stability the eigenfrequenices are {VLl^ + Vtl + 2Kaja^a + '^i^bjS^by^'^ 
and zero. As one of these frequencies, UJ2 = 0, is markedly different from those in the 
interaction free case, the effect of squeezing in the position space and the momentum space 
is much pronounced in the vicinity of the critical point. 
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To further elaborate this issue, we show the minimal value of during the evolution of the 
coupled condensates and the corresponding EOF as functions of the interspecies interaction 
strength k in Fig |21 which explicitly confirms that the degree of entanglement increases 
drastically as the interspecies interaction k is close to the stability limit given by (jl9|) . In 
fact, both quantities change noticeably once k, > k^- It is noteworthy that similar increase in 
entanglement has previously been found in quantum phase transition of a spin chain model 

Q 

VI. DISCUSSION 

In the present paper the entanglement between the two components of a BEG condensate 
trapped in a double- well is studied analytically in the low excitation limit with the HPT, and 
its accuracy is confirmed by comparison with the exact numerical solution. As demonstrated 
in previous sections, the degree of entanglement, gauged by the two- mode entanglement 
parameter, depends strongly on the interspecies interaction that can be varied with the 



current technology 
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19|. We expect that our work can be apphed to study entanglement 
in two-component condensates such as Rb-Rb and Rb-K mixtures. Specifically, our result 
shows that the two components of the condensate can remain in the tunneling phase and 
yet get strongly entangled as long as > k > Kg- 

To gain more physical insight from our result, we note that Jax and Jay respectively 
represent the population difference (measured by the operator 6ha = d^dij — dJdL) and the 
phase difference (measured by the relative phase operator 6(j)a) of the a-species condensate 
in the two wells, where a = a,b |3|]. Therefore, the two-mode entanglement parameter, 
C,t = n — c, where 

n = 2{(5n^(t))(50^(t))-[Re(5n,(t)50,(t))]2}'^' 

= 2 {{5nlit)){5^lit)) - [Re{SMt)SUt))?f' , (38) 
c = 2 [Re{{6na{t)6Ut)){Snb{t)SMt)) - (5n„(t)5n,(t))(50„(t)50,(t)))] , (39) 

indeed measure the correlation of 6na and 6(j)a- If the two components are entangled and 
therefore < 1, the fluctuations in the population difference and phase difference of the 
composite system are squeezed accordingly. 

On the other hand, in real experiments the temperature of the condensate is not exactly 
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zero |33|. Therefore, it is worth studying how the effect of finite temperature might affect 
the entanglement that could be built up in the system during its evolution when the initial 
state is a mixed state. Specifically, we consider an initial state that can be written as a 
product of two thermally equilibrium states maintained at a common temperature T, and 
the density matrix at t = is given by: 



p(0) = Pa® Pb, 



(40) 



where ioi a = a,b, 



Pa 



[1 - exp{-Qa/kBT)]exp{-Qadlaa/kBT). 



(41) 



It is a mixed stated and its two components are obviously separable at t = 0. The covariance 
matrix at t = 0, M(0), for this initial state is: 

/ 



M(0) 



2na + l 

2fi„ + 1 

2nfc + l 

2nfc + l 



\ 



where, as usual, the mean excitation number Ha is [exp{Qa/kBT) — 1] ^. If the initial 
temperature and the tunneling frequency are of order 10~^ K and 1 kHz respectively, which 



IS 



the mean excitation number may 



are typical values in current experiment situations 
reach order unity and can give rise to non-negligible effect on the entanglement parameter. 

It is well known that the Wigner characteristic function of a harmonic oscillator in a 
thennal ..xed .ate ,s stiH a Gaus.an f..t,o„ Q. So. ,t . .eadUy sKown tKat the 
system considered in our paper remains in a Gaussian state and one can use the two-mode 
entanglement parameter to study the entanglement between the condensates. Following 
the argument outlined previously, one can show that the covariance matrix M{t) at t > 
is given by U{t)M{0)U'^ [t), from which the two-mode entanglement parameter C,t can be 
obtained accordingly. In particular, for a symmetric two-component BEG with Qa = ^b, 
the covariance matrix M(t) is just U(t)U'^(t) multiplied by 2na + 1 (or equivalently 2nf, -|- 1) 
and the system is again symmetric. The two-mode entanglement parameter can then be 
used to determine whether the system is entangled or not and to evaluate the entanglement 
of formation EF{C,t) as well. In Fig. |3]we show the two- mode entanglement parameter C,t 
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and the entanglement of formation Ep as functions of time for a symmetric two-component 
BEC with the mean excitation number = 0.5. It is clearly manifested that a substantial 
amount of entanglement can still be achievable when k is sufficiently strong. Therefore, the 
existence of finite thermal effects does not readily preclude the occurrence of entanglement. 

On the other hand, for an asymmetric two-component BEC with Qa 7^ ^b, the covari- 
ance matrix M{t) is no longer symmetric with respect to the two interacting components. 
Therefore, despite that one can still make use of the parameter C,t to determine whether the 
system is entangled or not, it is no longer possible to apply the method developed here to 
obtain Ep even if < 1 and other numerical schemes have to be sought [s^]. 

In summary, although there might be some complications in quantifying the degree of en- 
tanglement if the initial state of the condensate is a mixed state, the two-mode entanglement 
parameter studied here can still give a necessary and sufficient condition of the separability 
of the two condensates. 
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APPENDIX A: EXPLICIT FORM OF THE VARIANCES 

The explicit expressions of the variances are given here for reference: 



{€) = E cos2cu,t + G, cos[a;i + (-l)'+^cu2]t + Hi} (Al) 



4 = 1,2 



= -FtT. {^^l^aF^ COS 2ujit - Q.G, cos[a;i + (-l)'+^^2]t + f^l^aHi} (A2) 



'a i=l,2 



iPl) = E {oo^F,cos2uj,t + {-iy+^UJiUJ2GiCOs[uJi + {~iy+^i02]t 




(Pi) = 777^ E {f^lujfQ,FiCOs2uJit + {-iyuj,uj2QbG^cos[uj, + {-iy+^uj2]t 
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1 
2 



(A4) 



R-e((g„g6)) = {'2l^iFiCos2i^it+ + ii2)GiCos[i^i + {-iy+^uj2]t + 2iiiHi}, (A5) 

^ i=l,2 

R'e((paPb)) = ^ ^ {2/Xia;,^FjCos2wii 

+(-1)*+Vc^2(a*i + At2)G'i cos[u;i + (-l)*+^u;2]i + 2//^^;,^//^}, 



i=l,2 



+ 



Re((g6p„)) = — — {"^l^iUJiFi sin 2ujit 



20. 



a. i=l,2 



+ [/X2a;i + (-l)^+Vi^2]G',sin[a;i + {-iy+'u2]t}, 



where 



H, = 



4(/xi- 


/^2)2 


1 




2(//i- 




~9 

/T 




4(/^i- 





1 nl nl 

IH ^ — 



02 



a;ia;2 



a;ia;2 



, 1 fi« 
1 H \ -\ — 



(A6) 
(A7) 

(A8) 

(A9) 

(AlO) 

(All) 
(A12) 
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FIG. 1: The time evolution of (solid line, left scale) and the EOF Ep (dashed line, right scale) 
are shown for the symmetric case with Na = Nf, = 400; Ka = Kb = h = = 50; and (a) 
K = 0.50; (b) K = 1.00; (c) k = 1.12. To gauge the accuracy of the HPT, the exact numerical 
solution of is shown by the empty circles. 
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FIG. 2: The time evolution of (solid line, left scale) and the EOF Ep (dashed line, right scale) 
are shown for the asymmetric case with Na = Nj, = 400; Ka = 0.75, Kf, = 1; fia = 34.5, fib = 50; 
and (a) k = 0.50; (b) k = 0.875; (c) k = 0.965. To gauge the accuracy of the HPT, the exact 
numerical solution of is shown by the empty circles. 
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K 



FIG. 3: The maximal achievable values of (empty circles, left scale) and the EOF Ep (crosses, 
right scale) during time evolution are shown for different values of the interspecies interaction 
strength k. The condensate is an asymmetric one with Na = Nf, = 400; Ka = 0.75, Kb = 1; and (a) 
na = 34.5, = 50, Kc = 0.969; (b) = 172.5, = 250, Kc = 1.385. The value of Kc is shown 
by the vertical dashed line in each diagram. 
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FIG. 4: The time evolution of (solid line, left scale) and the EOF Ep (dashed line, right scale) 
for a symmetric thermal BEC with Na = Nj, = 400; = = 1; J7a = $7;, = 50; fia = fif, = 0.5; 
and (a) k = 0.50; (b) k = 1.00; (c) k = 1.12. 
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